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by nonvanishing nonlocalizable information under the zero-way CLOCC protocol are 
expensive in computational cost. We define and evaluate economical measures on the 
basis of a new class of maps, eigenvalue-preserving-but-not-completely-eigenvalue- 
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q-i completely-positive (PnCP) maps that have been commonly used in the entangle- 



preserving (EnCE) maps. The class is in analogy to the class of positive-but-not- 



ment theories. Linear and nonlinear EnCE maps are investigated. We also prove 
subadditivity of the measures derived from the maps. 
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I. INTRODUCTION 



Nonclassical correlation of a bipartite system is an essential resource to perform quantum 
information processing.-^ Entanglement, namely the degree of inseparability, of a system is 
the most well-known nonclassical correlation. Besides the entanglement paradigm, several 
different paradigms 3-5 have been proposed in which the set of the states with nonclassical 
correlation includes certain nonentangled states. Unlike entangled states that are defined 
on the basis of the impossibility of preparation under local operations and classical com- 
munications (LOCC),- nonclassically correlated states in the different paradigms have been 
defined on the basis of post-preparation stages. 

Historically, nonlocality about locally nonmeasurable separable states was discussed by 
Bennett et al.— Later, in the context of system-apparatus correlation, a measure called quan- 
tum discord was defined by Ollivier and Zurek, 4 which is a discrepancy of two expressions 
of a mutual information that are equivalent in the regime of classical information theory. 
Recently, the term quantum discord often indicates the minimized one over the possible local 
(orthogonal) projection sets and is widely used as a measure of nonclassical correlation. As 
a typical example that justifies the computational power of nonclassical correlation other 
than entanglement, it was reported by Datta et al. 7 that the Knill and Laflamme's fast esti- 
mation of a normalized trace of a unitary matrix, which uses a single pseudo-pure qubit and 
the remaining qubits in a maximally mixed state, 8 exhibits a large quantum discord and 
vanishing entanglement. Another well-known definition of nonclassical correlation is the 
quantum deficit proposed by Oppenheim et al.,— which is nonlocalizable information under 
the closed LOCC (CLOCC) protocol that allows only local unitary operations and opera- 
tions to send subsystems through a complete dephasing channel. Among different setups 
they considered,- the zero-way setting — namely a setting in which the players are allowed 
to communicate under CLOCC only after local complete dephasing possibly subsequent to 
local unitary operations — was connected to a mathematically simple classical/nonclassical 
separation. 

The quantum deficit for a density matrix p AB of a bipartite system AB is defined as^ 

mm Ae cLocc[SvN(Tr B Ap AB ) + 5 vN (Tr A Ap AB )] - S v ^(p AB ), 

where SVnG) is the von Neumann entropy. Here, the system may be kept by a single player 
after the process A, i.e., TreAp AB or Tr A Ap AB possibly becomes a null state. The entropy of 
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a null state is defined to be zero. In case of the zero-way setting, the minimum is obtained 
when a single player possesses the total system after the process. Thus the zero- way quantum 
deficit is equal to a minimum discrepancy between S v ^(p AB ) and <S , v n[(V a <8> V B )p AB ] where 
V A ' B is a dephasing operation (acting on a local system) deleting the off-diagonal elements 
of a target matrix on a certain local basis; the minimum is taken over all local bases. 

The zero-way quantum deficit vanishes if the state has a product eigenbasis. A quantum 
bipartite system is called (properly) classically correlated 9 if and only if it is described by a 
density matrix having a product eigenbasis (PE), 

d A ,d B 

p£i = £ e jk \v A )(vf\^\v B )(v B \, (1) 
j,k=i,i 

where d A (d B ) is the dimension of the Hilbert space of A (B), is the eigenvalue of p 
corresponding to the eigenvector \v A ) <g> \v B ). Thus, a quantum bipartite system consisting 
of subsystems A and B is nonclassically correlated if and only if it is described by a density 
matrix having no product eigenbasis. 

Using this simple classical/nonclassical separation, other measures^ 1 ^ were later pro- 
posed. In particular, Piani et al.— recently designed a measure which vanishes if and only if 
a state has a product eigenbasis. It is in a similar form as quantum discord^ and defined as 
a distance of two different quantum mutual informations that is minimized over local maps 
associated with local positive operator-valued measurements.— 

A pending problem is that the original nonlocalizable information and the Piani et a/.'s 
measure both require expensive computational tasks to take minimums over all possible local 
operations. A similar difficulty exists in Groisman et al's measure^ which is a discrepancy 
between an original state p AB and the state after dephasing under an eigenbasis of Trep AB <8> 
TrAp AB (called Schmidt basis). In fact, the obviously classically correlated state ( 1 00) (00 1 + 
|ll)(ll|)/2 is mapped to 1/4 by dephasing when the improper Schmidt basis {|0), |1)} <8> 
{|±)} with |±) = (|0) ± |l))/\/2 is chosen while it is mapped to itself when a proper 
Schmidt basis (the computational basis in this case) is chosen. Thus a minimization over 
possible Schmidt bases is required. The recently-proposed measurement-induced disturbance 
measure^^ is a variant of the Groisman et al.'s measure; the same problem exists. In 
general, the measures involving a minimization over local operations are intractable in view 
of computational cost. 
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In our previous contribution,— an entropic measure G based on a sort of game to find 
the eigenvalues of a reduced density matrix from the eigenvalues of an original density 
matrix was proposed, in the context of m-partite m-split nonclassical correlation. This 
measure can be computed within a finite time although it does not have a perfect detection 
range. Its computational cost is indeed less than those for the intractable measures, but still 
exponential in the dimension of the Hilbert space. 

Here we introduce a new class of maps to define measures with improved computa- 
tional cost, in the context of bipartite splitting. It is the class of eigenvalue-preserving- 
but-not-completely-eigenvalue-preserving (EnCE) maps. We find it analogous to the class 
of positive-but-not-completely-positive (PnCP) maps^~— that are popularly used for detec- 
tion and quantification of entanglement. The idea of introducing the class EnCE was briefly 
mentioned in our previous contribution.— Here, we give mathematically strict definitions 
and show the fact that any linear EnCE map is a concatenation of unitary and anti-unitary 
operations. Thus the restriction of the theory is clarified. We further introduce a measure 
using a nonlinear EnCE map in order for achieving a wider detection range. 

The measures we propose here on the basis of EnCE maps are not as strong as those 
using an infinite number of trials, in the detection range of nonclassically correlated states. 
We propose a simple way to relax this drawback: The detection range is improved by 
introducing an average of multiple measures whose detection ranges are mutually different. 
This approach is described in Sec. IIVDI 

This paper is organized as follows. We start with an overview of the conventional theory 
of PnCP maps in Sec. [Til We then define and evaluate new classes of maps in Sec. II III The 
measures are defined and their properties, such as subadditivity, are verified in Sec. IIVI First 
non-subadditive measures are introduced in Sec. IIV Al Second subadditive measures are 
introduced in Sec. IIV Bl with the proof of the subadditivity. The computational complexity 
of the subadditive measures is investigated in Sec. IIV CI A simple way to relax the drawback 
in the detection range of the measures is shown in Sec. IIVDI Section [V] summarizes our 
results with some remarks. 
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II. CONVENTIONAL THEORY OF PNCP MAPS 



Quantum physics is governed by completely positive (CP) maps. Any map which is 
not CP (nCP) is considered to be physically unfeasible. There is, however, a class of nCP 
maps which are useful for characterizing entanglement. These maps are in the class of 
positive-but-not-completely-positive (PnCP) maps. It has been more than a decade since 
the Peres- Horodecki criterion opened the mathematical study of PnCP maps.— _ — A PnCP 
map Ap n cp is positive when acting as a global operator but nonpositive when acting as 
I eg) A PnCP on a system. It maps a separable state p sep = ^ Wipf <S> pf of a bipartite 
system AB with nonnegative weights Wi to a certain (physically feasible) state, while it does 
not necessarily map an inseparable state to a positive Hermitian matrix. Thus one finds a 
density matrix inseparable if one detects a negative eigenvalue of the matrix obtained after 
applying / Cg) Ap n cp to the density matrix. The PnCP map theory has gathered a broad 
interest in relation to detecting entanglement (See, e.g., Ref . fl) . 

One might be curious to find an analogue of the PnCP map theory to detect nonclassical 
correlations often defined in different ways^- than that of entanglement. We pursue the 
analogous theory to detect nonclassical correlation in the context of classical/nonclassical 
separation given by the existence/absence of a product eigenbasis of a bipartite state. 

III. INTRODUCTION OF UNCONVENTIONAL CLASSES OF MAPS 
AND THEIR USE 

We aim to introduce an analogy of the PnCP map theory to the present paradigm of 
classical/nonclassical separation. For this purpose, we define our new classes of maps. Let 
us start with a linear map theory. 

Definition 1. An eigenvalue-preserving (EP) map Agp is a map acting on a general d x d 
density matrix p = Ylt=i e k\ v k){ v k\ (here, and l^)^] are the kth eigenvalue and the 
corresponding projector, respectively) such that 




k=l 



k=l 
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where {\vk)}k and {\v'k)}k are both complete orthonormal systems (CONSs). The dimension 
of p' is equal to d. 

Alternatively, we may define the EP map in the following way: 

Definition 2. An EP map Agp acting on a quantum system S is a bijection between the set 
of projectors {\v k) (v k\}t=i generated from the vectors \vk) of a CONS to the set of projectors 
{\ v ' k) (v' k\}t=i generated from the vectors \v'k) of a CONS for any CONS {\vk)}f = i of the 
Hilbert space of S. 

A class of EP maps analogous to CP is defined as follows. 

Definition 3. An EP map A is a complete EP (CEP) map if and only if / <g> A is also an 
EP map for identity map I of arbitrary dimension. We denote such A as A C ep- 

Observation 1. One of the simplest CEP maps is U(d) : p — > upv) where u is an element 
of a flag manifold U(d) = U(<i)/U(l) xd (here, U(d) is the ci-dimensional unitary group). 

We now define the PnCP analogy in the following way. 

Definition 4. An EP map A is an EP-but-not-completely-EP (EnCE) map if and only if 
there exists an identity map I of some dimension, for which / <g> A is not an EP map. We 
denote such A as AehCe- 

We have defined a class of EnCE maps. As is analogous to the usage of a PnCP map, the 
usage of a linear EnCE map is to find a certain change of eigenvalues of a density matrix by 
applying the map to a local subsystem. This is based on the following proposition. 

Proposition 1. Both I A <g> A| nCE and A EnCE <8> I B preserve the eigenvalues of a density 
matrix of a system AB if the density matrix has a product eigenbasis. 

Proof. Let the density matrix with a product eigenbasis {1^)1^ )}j k =\ i be 

d A ,d B 

^ AB = E ^><^k B ><^ B i, 
j,k=i,i 

where is the (jk)ih eigenvalue corresponding to the (jTc)th projector \vf)(v^\ ® \vf)(v®\. 
By the definition, any linear EP map acting as a local operation should map an eigenbasis of 
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the reduced density matrix of a target subsystem to another CONS. Therefore, it is obvious 
that 

d A .d B 

(I A ® A| nCE )a AB = e jk \vf)(vf\ ® \v' B k )(v' B k \, 

j,k=l,l 

where {|?/ B )}fc is a CONS of the Hilbert space of B, which may be different from {|f B )}fc. 
It is trivial to show the same proof applies to A EnCE £g> I B . □ 

Corollary 1. A density matrix p AB has no product eigenbasis if either (I A <S> A EnCE )p AB or 
(A EnCE <S> I B )p AB has eigenvalues different from those of p AB . 

Proof. This is the contraposition of Proposition HJ □ 

There is, however, a restriction in the type of the linear EnCE maps according to the 
following proposition. This restriction is later relaxed by a nonlinear EnCE map. 

Proposition 2. Any linear EP map can be decomposed into unitary transformations and a 
transposition. Hence any linear EnCE map can be decomposed into unitary transformations 
and a transposition. 

Proof. Consider a linear EP map Aij n and two pure states \x) (x\ and \y) (y\ (\x) and \y) can be 
nonorthogonal to each other) . Let \x')(x'\ = Ai in (|x)(x|) and \y')(y'\ = An n (\y) (y\) . Consider 
a state r = \x)(x\ + \y)(y\ represented under a certain CONS. As A lin changes this CONS 
to a certain CONS, (x\r\x) is equal to (x'\An n (T)\x'). This suggests that = | {a;' | j/') | 

since Ai in is a linear map. Note that this is true for any linear EP map Ai in and any two pure 
states \x) and \y). Therefore, by Wigner's unitary-antiunitary theorem,— there are only two 
possible types for Ai in , namely, unitary and antiunitary^ transformations acting on a target 
density matrix. Hence the proposition holds. □ 

Observation 2. A unique nontrivial linear EnCE map is the transposition At, according 
to the above proposition. As an example of detecting a nonclassical correlation, consider 
the density matrix of a two-qubit pseudo-entangled (PS) state, 

p ps = (l-p)//4 + p|V)(V| (2) 

with \vp) = (| 00) + 1 1 1) )/ -\/2 and < p < 1. It has a nondegenerate eigenvalue (l + 3p)/4 and 
a degenerate eigenvalue (1 — p)/4 with multiplicity 3. Its partial transposition, (7® A T )p ps , 
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has a nondegenerate eigenvalue (1 — 3p)/4 and a degenerate eigenvalue (1 + p)/4 with 
multiplicity 3. These two sets of eigenvalues are different for p > 0, indicating the existence 
of a nonclassical correlation. 

It should be noted that having different eigenvalues after partial transposition is a suffi- 
cient but not necessary condition for a state to have no product eigenbasis. For example, a 
2-qubit state 

Po+ = i(|00)(00| + | + +>(+ + |) (3) 

with |+) = (|0) + |l))/\/2 has no product eigenbasis because |0)(0| and |+)(+| cannot be 
diagonalized simultaneously. It is clear that the partial transposition does not change the 
state and hence it does not detect a nonclassical correlation. We may use, instead, the 
nonlinear map V x defined later in order to detect a nonclassical correlation of this state. 

The use of a linear EnCE map for detecting nonclassical correlation is intuitive and 
computationally easy as we have seen. There is, however, a case where the limitation of the 
linear EnCE maps is clearly observed as described below. 

Remark 1. There are states called one-way classically correlated (1WCC) states,- in the 
form 

PlWCC = l iX )^ X ! ®°i 
i 

with \i x ) a CONS of x = A or B and y the remaining system; of (unnormalized) density 
operators acting on y, dependent on the index i. Such a state may have no product eigenbasis 
but testing a change in the eigenvalues under J x £g> Ag nCE for a single side is not enough to 
detect it. Therefore we need to test for both (x, y) = (A, B) and (x, y) = (B, A). 

Proposition 3. One cannot detect a nonclassical correlation of a one-way classically cor- 
related (1WCC) state using a linear EnCE map. 

Proof. It is easy to find that applying a partial transposition to piwcc results in either 
{U*®P) Piwcc (U*®Fy or (U*®Fy p* lwcc {U*®P) with £/„ = In addition, 

any partial unitary transformation preserves the eigenvalues of piwco By Proposition [2, 
the proof is completed. □ 

This proposition suggests that we need to search for nonlinear EnCE maps for a wider 
range of detection than that of linear ones. The definition involving both linear and nonlinear 
ones should be newly given in consistent with Proposition [TJ 
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Definition 5. An EnCE map AehCE (that can be nonlinear) should have the following 
properties. 

(i) For any density matrix p = J2k=i e k\ v k)( v k\, 

A E nCE : P J2k=i e k\ v k)( v k'\ where {\v k )} and {\v k ')} are both CONSs; e k are the 
eigenvalues. 23 

(ii) For any bipartite density matrix ppf with a product eigenbasis, written as (PQ), 

I A <g) A B nCE : ppf i — y p AB where p AB is a square matrix with the set of the eigenvalues 
same as that of p AB . 

(iii) For some bipartite density matrix <r AB having no product eigenbasis, -f A ®A B nCE maps 
it to a square matrix with the set of the eigenvalues different from that of a AB . 

We find that there is, in fact, a useful nonlinear EnCE map. To define it, we first introduce 
the specially-designed nonlinear map T x . 

Definition 6. A nonlinear map r x acting on a (possibly unnormalized) quantum state p is 
defined as follows. 

T x :p^ v /(pp- 1 )(h.c.) = p-, 



I A <8> T B : p AB ^ v/{p AB [/ A ® (Tr A p AB ) 1 - 1 ]}{h.c.}, 

where iGR; the square root is positive and h.c. stands for Hermitian conjugate (conjugate 
transpose). 

[Here, the inverse of a density matrix p = J2 k ek\vk)(v k \, with (e^, \vk)) the pair of an eigen- 
value and the (normalized) corresponding eigenvector, is defined as p _1 = J2k Cfc ^o c fe \ v k) ( v k\ 

This is an extension of the xth power of a matrix. Note that (I A £§> T B )p AB is a quan- 
tum state (positive Hermitian matrix) because, for positive Hermitian matrices A and B, 
(AB)(h.c.) = ABBA is a positive Hermitian matrix. 

A nonlinear EnCE map is now defined by using T x . 

Definition 7. A nonlinear EnCE map V x acting on a quantum state p is defined as follows. 

^:p AB ^r 1/a r> AB = p AB , 

i A ® p b : p AB ^ {i A ® r B j(/ A ® r B ) P AB , 

where igR,i^L 



Of course, (7 A <g> V B )p AB is a quantum state (the trace is not preserved in general). 
The map V x is useful for detecting nonclassical correlation because we have the following 
theorem. 

Theorem 1. The equations (I A ®V B )p AB = {V A ®I B )p AB = p AB hold if p AB has a product 
eigenbasis. 

Proof. For a bipartite state with a product eigenbasis (PE), p AB = c ij\ u i) A ( u i\ ® 

\ v j) B ( v j\i we nave Tr APpE = J2l(J2k C kl)\ v l) B ( v l\- ^US 

= (/ A r B )p AB = ^pABj/A^^AB),-!]}^ ^ 
= Ej, /(^o/0') X '~ 1 E i C ij |M i ) A (M i | ® |%} B (%| 

with /(j) = Cfcj. For this matrix, we have 

Tr A p^ = J2 f(ty- l f(t)\u t ) B (u t \= f(t) X \ut) B (u t \. 

t, f(t)ytO t, f(t)ytO 

Thus 

= Ei, /ovo/OT '' E, r v "/) A (», ® |^»,| 
= p ab . 

This proves that (I A <g> P B )ppf = p A |. It is easy to show that (P A g> J B )ppf = Pp| in the 
same way. □ 

Here is a simple example to use this map for detecting nonclassical correlation. For 
the bipartite state po+ which has been introduced in (EJ), (J £g> V2)po+ has the eigenvalues 
(approximately) 0.826, 0.375, and (strictly) (with multiplicity two) which are different 
from the eigenvalues of po+, 3/4, 1/4, and (with multiplicity two), except 0's. Therefore 
Po+ has no product eigenbasis. 

IV. QUANTIFICATION OF NONCLASSICAL CORRELATION 

For the next step, we define a measure of nonclassical correlation based on the theory 
of EnCE maps we have seen. Note that, according to the definition, the set of classically 
correlated states is a nonconvex subset of the set of separable states. Thus it is not motivating 
to impose convexity on a measure of nonclassical correlation. We may, however, impose a 
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family of additivity properties.— In particular, subadditivity is assessed in the following. We 
begin with non-subadditive measures and later introduce subadditive measures. A strategy 
to extend the detection ability of subadditive measures is described. 

A. Non-subadditive measures 

We first define a non-subadditive measure of nonclassical correlation for a given EnCE 
map AEnCE as follows. Suppose we want to quantify a nonclassical correlation of a bipartite 
system AB described by a density matrix p AB . Then, we may consider the quantity with 
subscript R (L) indicating that the right (left) component is acted by AehCe: 

£>R,L(A E nCE,p AB ) = ^ |e s - e' s |, 

s 

where e s 's are the eigenvalues of p AB while e' s 's are those of (J A ® A B nCE )p AB for "R" 
[(A EnCE ®/ B )p AB for "L"]; e s 's and e' s 's are aligned, say, in the descending order. In general 
e' s may possess an imaginary part while it does not when we use the transposition At or 
the map V x defined in Definition [7] for A En cE- It is obvious that D vanishes if p AB has a 
product eigenbasis. 

We can easily calculate -D R L (A EnCE , p AB ). For example, it is easy to calculate^ D R (A T , p ps ) = 
2p for the two-qubit state p ps defined in fl2]). 

Another simple example is the quantification of a nonclassical correlation for the bipartite 
state po+ which has been introduced in ([3]). The eigenvalues of po+ are 3/4, 1/4, and (with 
multiplicity two). The quantity _Dr(At,Po+) vanishes because (I £g> Ax)po+ = Po+- In 
contrast, D R (V2, Po+) does not vanish: as we have computed in an example in the previous 
section, (I ®V2)po+ has the eigenvalues (approximately) 0.826, 0.375, and (strictly) (with 
multiplicity two). Thus D R (V 2 , Po+) ^ 0.201. 

The last example is to clarify that D RL is not subadditive. For \ip) = ( 1 00) + |ll))/\/2, 
(J A <S> A B )|?/>) AB (-?/>| has the nondegenerate eigenvalue —1/2 and the degenerate eigenvalue 
1/2 with multiplicity three. The eigenvalues of are 1 and with multiplicity three. 

Thus -Dr(A t , \ip)(ip\) = 2. It is now easy to find the eigenvalues of (J AC <g> A^ D )|^) AB (^| ® 
\ijj) CB , which are —1/4 with multiplicity six and 1/4 with multiplicity ten. This results in 
D R (A T , ® = 9/2. This value is larger than 2 x D R (A T , \ip)(ip\) = 4, indicating 

that subadditivity does not hold. 
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B. Subadditive measures 



It has been shown that the measures -Dr,l introduced above are neither additive nor 
subadditive. Additive or subadditive measures are desirable if one needs to compare systems 
with different dimensions. Here, subadditive measures are introduced. Let us formally begin 
with the definition of subadditivity.— 

Definition 8. Let F(p AB ) A | B be a measure of correlation between subsystems A and B of 
a bipartite system AB, where A|B denotes splitting between A and B. Then, F(/? ab )a|b is 
called a subadditive measure if and only if the relation F(p AB <g) ct cd )ac|bd < -FXp AB )a|b + 
F(a CD ) C |D holds. 

We find that the following quantities Qr and Ql satisfy the subadditivity condition if we 
choose the map AehCE properly. We define them as 

Qr,l(AeiiCE, P AB ) = - log 2 V 7 ^) > ( 4 ) 

where e s 's are the eigenvalues of p AB and e^'s are the absolute values of the eigenvalues of 
(I A g> A| nCE )p AB for "R" [ (AgnCE ® I B )p AB for "L"]; e s 's and e s 's are both sorted, say, in 



descending order; N = v/J^s is a normalization factor which guarantees Qr,l > 0. The 
measures Qr,l vanish if {e s } = {e s /N 2 }. As for subadditivity, we can prove the following 
proposition. 

Proposition 4. The measure QR(AEnCE> P AB ) is subadditive if the set of the eigenvalues of 
(I AC (g) A|° CE )(p AB Cg> cr CD ) is given by {ajb k } jk with the eigenvalues aj of (I A <g) A EnCE )p AB 
and the eigenvalues b k of (I c ® A EnCE )cr CD . 

Proof. The proof consists of two steps (i) and (ii). 

(i) Consider the two sequences {pi}f =1 and {qi}f =1 of nonnegative real numbers pi and g;. 
Suppose they are sorted: pi > p 2 > ■ ■ ■ > Pd and q 1 > q 2 > ■ ■ ■ > q<i- Then, the fidelity 
Y2t=i y/PiQi f° r these sorted sequences is larger than that for any two unsorted sequences 
whose entries area's and q^s, respectively. This is because, for real numbers a±, a 2 , bi, and 
b 2 such that a x > a 2 and b\ > b 2 , the relation a\b\ + a 2 b 2 > a\b 2 + a 2 bi holds. 

(ii) Let us write the eigenvalues of p AB as <%,• and those of cr CD as b k . The fidelity F' = 
Sj=i ^2k=i \j ( a jbk)(ajb k )/N with N = \ f^~ h Sjb k involves two possibly unsorted sequences 
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{{ a jbk)} and {(aj&fc)}; these are unsorted in general even when {aj}, {bk}, {Sj}, and {b k } 
are individually sorted. Let us write the fidelity after sorting {(oj&fc)} and as F. 

Then, F' < F < 1 holds according to the fact (i). Therefore, < — log 2 F < —\og 2 F' = 



- log 2 F a -log 2 F b holds with F a = J2j \/ a j a j/ Jj2j a j and F b = Y.k\ b kh/yJ2k b k, where 



It is trivial to find a similar condition for Ql(AeiiCe, p AB ) to be subadditive. In addi- 
tion, it is clear that Qr,l vanish if p AB has a product eigenbasis. These measures are a 
sort of logarithmic fidelity and are reminiscent of logarithmic negativity.—"— We find that 
choosing AehCE from the maps At and V x introduced in the previous section satisfies the 
condition of Proposition H] as we prove below. As for other additivity properties, Qr,l is 
not additive or weakly additive in general owing to sortings of the eigenvalues. This is 
clear from the following example: For the state p ps defined in (T5]), with p set to 1/3, we 



have Qr(A t ,P AB =1/3 ® p£ s b =1/3 )ac|bd = -log 2 (5/18 + v^/3) ^ 0.226; this is less than 
2 x Qr(A t ,p ab =1/3 ) A |b = -21og 2 (V6/6 + V2/3) ~ 0.370. 



As we have mentioned above, one map that makes the measure Qr,l subadditive is the 
transposition At- The subadditivity is easily verified according to the fact that J AC £g> 
A BD = (J A £g> A B )(/ C eg) A°). The condition on the set of eigenvalues stated in Proposition 
H] is obviously satisfied. In addition, the measures Qr,l(At, p AB ) are invariant under local 
unitary operations. Its invariance under local unitary operations (say, U B ) follows from 



(J A ® A B )(/ A ® U B p AB I A ® W B ) = (I A ® U* B )(I A ® A B p AB )(I A ® U** ). 

It is also easy to find that V x , the map defined in Definition [TJ makes Qr,l subadditive. 
Because TtacP AB ® P CD = (TrAp AB ) <S> (Trcp CD ), we have 
(J AC <g> r BD )(p AB <g> a CD ) = v/{p AB [/ A <8> (TrAP^) 1 '- 1 ] ® a CD [J c ® (Tr c a CD )-- 1 ]}{h.c.} 



Hence, (7 AC g> P BD )(p AB ® (x CD ) = (7 A ® P B )p AB ® (7 C <g> P B )a CD holds. Therefore, by 
Proposition HI we find that Qr(V x , P AB ) is subadditive [and Ql(P x , p AB ) either]. In addition, 
it is invariant under local unitary operations because (J A ' B <g> r B,A )p AB does not depend on 
the choice of a local basis, by its definition. 

The problem is that Qr and Ql are different in general. To solve this problem, we suggest 
using the average 





{aj}, {bk}, {aj}, and {bk} are individually sorted. 



□ 



(/ A ® r B )p AB ® (/ c ® r B )a CD . 



Q(A En cE,p AB ) 



Qr(AeiiCE7 P AB ) + Q L (AEnCE, p AB ) 

2 



(5) 
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This becomes subadditive and invariant under local unitary operations if both Qr and Ql 
are subadditive and invariant under local unitary operations. It is easy to find that At and 
V x are both useful for this purpose. 

C. Computational complexity 

An advantage of using the measure Qr,l defined in (TjJ and the measure Q defined in (jSJ) 
is their relatively small computational cost when A T or V x is chosen for A EnCE . 

Consider a bipartite system with the dimension d A (d B ) of the Hilbert space of its sub- 
system A (B). For a density matrix p AB , (I (g) At)p ab is computed with 0(d A2 d B2 ) basic 
floating-point operations. This is less than the cost of diagonalization of (leg) A^)p AB , which 
takes 0(d A3 d B3 ) basic floating-point operations. The computation of (/ ®V x )p AB is a little 
expensive because it involves a square root of a matrix. The cost of computing (I <g> V x )p AB 
is 0(d A3 d B3 ) : same as the cost of diagonalizing (/ £g> V x )p . 

Once the eigenvalues of (7® A Eii ce)p AB is computed, it takes only 0(d A d B ) basic floating- 
point operations to compute Qr,l and Q. Therefore, the time complexity of computing these 
measures is 0(d A3 d B3 ) when At or V x is chosen. 

D. Extending the detection range 

One might be curious if V x , the map defined in Definition [7J is more useful than the 
transposition At in detecting nonclassical correlation by using the measure Q defined in ([5]) . 
First, Q(A T ,p AB ) vanishes while Q(V 2 ,Po B ) ^ 7.00 x 10" 3 for the state p + defined in 
Second, Q(At, \ip} AB = 1 while Q(V x ,\ip} AB (ip\) vanishes for |^) = (|00) + |11))/V2. 
Therefore, generally speaking, Q(V X , p AB ) is neither stronger nor weaker than Q(At,p ab ). 
One may further claim that V x is not very useful because it vanishes for the Bell state. 
Nevertheless, this is not a serious drawback as we have a quick solution as follows. 

There is a way to utilize these measures to produce a stronger measure. Suppose we 
have non-negative, subadditive, and local-unitary-invariant measures Mi, ...,Mn- Then, the 
weighted average w^M^ with Wk > is also a measure which is non-negative, subadditive, 
and invariant under local unitary operations. It detects nonclassical correlation for the states 
for which any one of Mi, M^ is nonvanishing. 
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Thus we easily produce the stronger measure 

w T Q(A T ,p AB ) + ^w k Q(V Xk ,p AB ) 
k 

with Xk G R, Xfc 7^ 1, and tOx, Wk > 0. This measure does not vanish for p Q+ and |^)(^|. 

V. CONCLUDING REMARKS 

We have seen several different usages of the EnCE map theory. We believe that this 
theory works as a useful template to detect and quantify nonclassical correlation based on 
the Oppenheim-Horodecki separation of classical/nonclassical correlations. The EnCE map 
theory has been constructed in analogy to the PnCP map theory in the present paper. 
One important difference between these theories is that the class of EnCE maps includes 
nonlinear EnCE maps. This is because linear EnCE maps are very limited due to the fact 
that any linear EP map can be decomposed into unitary operations and a transposition 
(Proposition [2]). Nonlinearity of a map is not a significant problem as far as I <g> A En c E and 
A-EnCE ® I are defined appropriately for an EnCE map A EnCE in the way that / <g> A EnCE and 
AeiiCE <E> I preserve the eigenvalues of any state that has a product eigenbasis. 

On the basis of the EnCE map theory, we have defined two subadditive measures, 
Q(Ax,p AB ) and Q(V X , p AB )- These are neither stronger nor weaker to each other in the 
detection range, and not so strong as the measure by Piani et al.— that is perfect in the 
detection range albeit intractable in computational cost. The advantage of our measures is 
the complexity: they are calculated within polynomial time in the dimension of the Hilbert 
space. We have shown a way to relax the drawback of the detection range; their weighted 
average is stronger than themselves and remains subadditive as shown in Sec. IIVDI A 
certain optimization over the weights and the choices of x's will be investigated in future 
work. 

One might be curious about an extension of the measures to multipartite splitting. This 
is achieved by taking a minimum, maximum, or average of a measure over all possible 
bipartite splittings of the multipartite system. In considering the possible combinations 
of subsystems for a bipartite splitting, we should be careful about the fact that having 
product eigenbases for A|BC splitting and AB|C splitting does not imply having a product 
eigenbasis for A|B|C splitting. A typical example is the state (|000)(000| + |1 + + l|)/2 
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with |+) = (|0) + |l))/v2. This state does not have a product eigenbasis for A|B|C splitting 
while it has for A|BC and AB|C splittings. A proper claim is that having product eigenbases 
for the A|BC, AB|C, and AC|B splittings implies having a product eigenbasis for A|B|C 
splitting. More generally, an m-partite state p 1 — m has a fully product eigenbasis if and only 
if p l - m has a product eigenbasis for every possible bipartite splitting separating {1, ...,m} 
into two sets. The proof is given in Appendix |A] 

In summary, a comprehensive framework, called the EnCE map theory, to detect and 
quantify nonclassical correlation of a bipartite system has been proposed. The average loga- 
rithmic fidelity Q(AEnCE, p AB ) has been introduced as a subadditive measure for a properly- 
chosen EnCE map Aedce- It is computable within polynomial time in the dimension of the 
Hilbert space. A simple way to extend the detection range by a collection of measures has 
been developed. 
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Appendix A: Theorem on multipartite product eigenbasis 

Theorem 2. An m-partite state p l --- m has a fully product eigenbasis if and only if p 1 — m has 
a product eigenbasis for every possible bipartite splitting separating {1, ...,m} into two sets. 

Proof. It is trivial that p l - m has a product eigenbasis for every possible bipartite splitting 
if it has a fully product eigenbasis. 

Now we prove the converse. By lemma [1] introduced below, the density matrix has 
a product eigenbasis for the l|2|34...m splitting and that for the 12|3|4...m splitting. The 
latter fact implies that p 1 ■ m 's eigenbasis is a product of the eigenbasis of the reduced density 
matrix p 12 , that of p 3 , and that of p 4 - m . The former fact implies that the reduced density 
matrix p 12 has a product eigenbasis. Therefore, p l - m has a product eigenbasis for the 
l|2|3|4...m splitting. 

16 



Next, we use the fact that p x - m has a product eigenbasis for the 123|4|56...m splitting by 
lemmaHJ Now it is found that p l - m has a product eigenbasis for the l|2|3|4|56...m splitting. 
Using the same logic continuously, the converse is proved. □ 

Lemma 1. A tripartite density matrix p ABC has a tripartite product eigenbasis if and only 
if it has a bipartite product eigenbasis for each of all the bipartite splittings. 

Proof. It is trivial that p ABC has a product eigenbasis for any bipartite splitting if it has a 
tripartite product eigenbasis. 

Now we prove the converse. Having a bipartite product eigenbasis for any bipartite 
splitting implies that 

p ABC = Ei 3 a ij\ r i) A ( r i\ ® h') BC (*jl 
= hi\t k ) AB (t k \ ® h) c («/| 

= Y.mn C mn\Vm) AC (v m \ ® | W n ) B (w n | , 

where |rj) A , |s-,) BC , |tfc) AB , \ui) c , |f m ) AC , and \w n ) B are eigenvectors of the reduced density 
matrices of the indicated subsystems; a^-, b k[ , and c mn are eigenvalues of p ABC . 
This leads to that 

(i) An eigenbasis of p ABC is a product of an eigenbasis of TrBcP ABC and that of TrAp ABC . 

(ii) Matrix TrAp ABC is represented as 

Tr A p ABC ®l«l>°<ti|| 



with a B = Tr A \t k ) AB (t k \ and o' G m = Ti A \v m ) AC (v m \. 
From (ii), we find that 

Tr AP ABC K) B K) c = (E k b ky a B \w x } B )\u y } c 

= l^) B (Em C m^ (7 ml%) C )- 

This implies that (%2 k bk y o~ B \w x ) B ) = \w x ) B . Hence {\w x ) B \u y ) c } is an eigenbasis of TrAp ABC . 
This fact and (i) complete the proof. □ 
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